Entanglement is not very useful for estimating multiple phases 



O ■ 

o : 

s ■ 

<: 
in : 

(N . 
CM ' 

> : 
o . 

On ■ 

m 
o 

o 

22 

Oh! 
-I— > 



Manuel A. BallesteiQ 

Department of Mathematics, University of Utrecht, Box 80010, 3508 TA Utrecht, The Netherlands 

The problem of the estimation of multiple phases (or of commuting unitaries) is considered. This 
is a sub-model of the estimation of a completely unknown unitary operation where it has been 
shown in recent works that there are considerable improvements by using entangled input states 
and entangled measurements. Here it is shown that when estimating commuting unitaries, there is 
practically no advantage in using entangled input states or entangled measurements. 



I. INTRODUCTION 

A unitary operation is a map that transforms a density 
operator p on C d to another density operator p = UpoW 
on C d , where U E SU(d) is a d x d special unitary matrix. 
Suppose one is given a device that performs an unknown 
U. One can learn something about U by learning about 
how it transforms a known state po- In order to com- 
pletely determine a unitary operation one would need to 
know how it transforms a basis of C d plus some linear 
combinations thereof. This is known as quantum process 
tomography pj. More precisely, to estimate U one pre- 
pares many copies of the necessary input states and per- 
forms a measurement on the output that they produce. 
As a result, some classical data are obtained and from 
that one can estimate U. This is shown schematically in 
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FIG. 1: Quantum process tomography. 

Another approach (used in Refs. 0, IS EL 01) is to pre- 
pare a bipartite entangled input state po on C d £g> C d and 
then use one of the parts as input for U, while noth- 
ing is done to the other part, as shown in Fig. [21 The 
effect of the operation is to transform the state p$ to 
(U ® l)po(f^' ® !)■ This output state is then measured 
and estimated, and, since in this case there is a one-to- 
one relation between the output state and U, one gets an 
estimate for U as well. The advantages of this method 
with respect to quantum process tomography are that 
only one input state is needed, and that there is poten- 
tially a better accuracy in the estimation (if nonseparable 
measurements are used) 

In this paper a relatively less difficult problem will be 
studied, the estimation of unitary operations that com- 
mute with one another, that is, only a maximal Abelian 
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FIG. 2: Entanglement is used. 



C d ® C d model. 



subgroup of SU(d) is considered instead of the whole 
group. In this case, the number of unknown parameters 
decreases from d 2 — I to d — 1. This problem has already 
been addressed in Ref. @ where it is given the name of 
"multiple phase estimation" (MPE). One would like to 
know whether it is also advantageous to use an entan- 
gled input in MPE. In what follows the MPE model that 
uses entanglement (Fig. |5J will be referred to as MPEE 
and the one that does not use it (Fig. ^| will be referred 
to as MPEU. 

There are two things that need to be optimized here, 
the input state and the measurement that is to be per- 
formed. Therefore one needs a quantitative measure of 
how good an input state is and of how good a measure- 
ment is. 



II. 



PRELIMINARIES 



In this section, the necessary concepts of quantum 
Fisher information (QFI) and Fisher information (FI) 
will be introduced, and the quantum Cramer-Rao bound 
(QCRB) will be stated. The QFI and the FI will be used 
as measures of the performance of an input state and a 
measurement, respectively. The QCRB relates these two 
quantities in a nice way. 



A. QFI 

Suppose that the quantum state density matrix p on 
C d is parametrized by 9 € O C R p where p is the number 
of parameters. In our case p would be the output state 
and p = d — 1. Define the symmetric logarithmic deriva- 



tives (SLDs) Ai, 
satisfy 



, X p as the self-adjoint operators that 
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p,i{6) = de iP (e) = ±[p(0)Ai(0) + \i(6)p{6)]. 
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For pure states, p = \ip)(ip\, they simply are = 2pj. 
The QFI is defined as the p x p matrix with elements 

Huiff) = Retr b(0)Ai(0)Aj(0)] 

which for pure states reduces to 

77^(0) =Re<Z i (0)|Z J (0)) 

where |Z»(0)) = Xi(9)\ip(9)) . 

The |Z) vectors have a simple geometric interpretation. 
Suppose one has a pure state model parametrized by 
9 € C R. For simplicity, the parameter has been taken 
to be one dimensional. Denote the set of all state vectors 
by H = {\ip(9))\9 € 0} and the set of all density opera- 
tors by M = {p{9)\9 G 0}. State vectors and density op- 
erators are related by the map 7r : | ip) i— ► \tp) , H — ► .M . 
This map is many to one since 7r(|V')) = """(e^lV')) where 
<j) is an arbitrary real phase. This means that 7r -1 (/?o) is 
a circle in Ti and that a curve in M. is mapped through 
tt" 1 to a tube in Tt. Conversely, all curves that lie on the 
surface of the tube, are mapped through tt to the same 
curve in M.. Then out of all curves \ip{6)) € TL satisfying 
\ip{0)) — \ipo) and tt(\^(9))) = p{6) (p(9) is a given curve 
in M) there is a minimal curve \ip(9)), defined as the 
one that at every 9 goes from the circle tt" 1 [p(9)] to the 
circle ir~ 1 [p(9 + 59)] using the shortest path. \l(0)) is a 
vector pointing in the direction of that shortest path. It 
can then be calculated as \l(6)) = 2dg\tp(9)). 



B. FI 

Take a positive operator values measure (POVM) with 
elements Mi, . . . , M n . This POVM induces a probability 
distribution given by p^(9) — tr p(9)M^ the probability 
to obtain outcome £ if the parameter has the value 9. 
The Fisher information for this measurement is defined 
as the p x p matrix with elements 

Iij(M,e) = EmA^ \n Pi {9)d B] \n Pi {9)\. 

For an estimator 9 and a measurement M, locally unbi- 
ased at 9o 1 the (classical) Cramer- Rao bound is satisfied 

V(M,9 ,§)>I(M,9 )-\ 

i.e., the FI is the smallest variance that a locally unbiased 
estimator based on this measurement can have. This also 
means that, if one of the eigenvalues of 7 is zero, then the 
variance of the function of the parameters corresponding 
to that eigenvalue is infinity and therefore cannot be es- 
timated. 



1 This means that the expectation of the estimator satisfies 
EM,e o 0i) = &0i and d 0j E M ,80i)\e=0 o = 5y- 



If one has N copies of the quantum state and performs 
the same measurement on each of the copies then the 
FI of the N copies, 7^, satisfies I N (M,9) = NI(M,6) 
where I(M, 8) is the FI of one system. It follows that 

V N (M, 9 , 9) > I N (M, 9q)~ 1 = 7(M,0 o r7^- 

It is a well known fact in mathematical statistics that 
the maximum likelihood estimator (MLE) in the limit 
of large iV is asymptotically unbiased and saturates the 
classical Cramer-Rao bound. Moreover no other reason- 
able estimator (unbiased or not) can do better. 

C. QCRB 

The QCRB states that for any measurement M 

I{M,9)<H{9). (1) 

In other words, H(9) — 7(M, 9) is a positive semidcfinitc 
matrix. 

This bound is not achievable in general. A theorem due 
to Matsumoto states that for pure states the bound 
is achievable at 9 — 9q if and only if 

lm(h(9 a )\l 3 (9 Q )) =0. (2) 

If a model satisfies the above condition, it is said to be 
quasiclassical at 9q. Furthermore, if condition J5J holds, 
there is a measurement with p + 2 elements that achieves 
the bound. In fact, any measurement of the type 

M a = \b a )(b a \, a=l,...,p+l, 

m+l 

M p+2 =t-J2M a , 

£J (3) 
\b a ) = 2_^o af3 \mp), 

\ m k) = ^(i7~5) fc; |; ; ) ; \m p+ i) = 
l 

with o a (p+1) x (p+1) real orthogonal matrix satisfying 
^ achieves I(M, 9 ) = H(9 ). 
Now one can see why these quantities are a good mea- 
sure of the performance of input states and measure- 
ments. From V > I^ 1 /N > H^ 1 /N one can see that 
a good input state is one that achieves an H as large as 
possible and a good measurement is one that achieves an 
7 as large as possible. Since I and H are matrices the 
best input state and best measurement cannot always be 
decided unambiguously. This ambiguity will not be com- 
pletely resolved in the case of choosing an input state. 
However, in the case of choosing an optimal measure- 
ment there will be no ambiguity. In the next section, it 
will be shown that for every input state there is a mea- 
surement on the output state that achieves equality in 
Eq. Q). Therefore, an optimal measurement is one that 
achieves equality in the QCRB. 
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III. MPE IS A QUASICLASSICAL MODEL 

It will be shown here that both MPEE and MPEU are 
quasiclassical everywhere (for all 9 £ and for any 

input state. In particular, this means that only one input 
state is necessary in MPEU also. In this respect MPE is 
quite different from the estimation of a completely arbi- 
trary U. Actually, since MPEU can be considered as a 
special case of MPEE, one needs to show quasiclassicality 
only for MPEE. 

The MPEE model is 



p{9) = UgpoUl 



(4) 



where 



Furthermore, there are situations in which one has two 
input states p\ and p2 with QFI H± and H2, respectively, 
which satisfy neither Hi < if 2 nor Hi > H%. This is re- 
solved by maximizing a quantity like Tr GH, with G a 
real positive semidefinite matrix. By doing this one as- 
signs relative weights to the mean square errors of the 
different parameters. Furthermore, achieving maximum 
Tr GH ensures that no other input state can have a larger 
QFI. In this paper a particular choice is made: all param- 
eters are given the same importance, i.e., the input state 
will be chosen so that it maximizes Tr H. With this par- 
ticular choice it is possible to obtain nice analytic results 
and, since the trace of the QFI is parametrization invari- 
ant, the results obtained will not depend on the chosen 
parametrization. 



U e = exp ( i mT m I ® 1, 



MPEE 



Po on C d ® C d is a pure state density matrix, 2 and 
Ti, . . . ,T<j-i are selfadjoint traceless matrices that com- 
mute with one another. They are chosen so that they 
satisfy an orthonormality condition 



tr T m T n — S nm • 



(5) 



The SLDs are 



\ m = 2d m p{6) = 2i[T m ®±,p{6)] 



and 



trpA ro A n = 4{trp (T m T n (g>iytr[po(T m <S>l)po(T n <S>l)}}. 

It is easy to see that due to the commutativity of the T's 
this quantity is real and therefore the model is quasiclas- 
sical, i.e., it satisfies the condition J5J. Therefore for ev- 
ery 9 there exists a measurement M (which may depend 
on 9) such that I(M,9) = H(9). Furthermore, if one 
has a large number N of copies, performs the same (op- 
timal) measurement on all N, and calculates the MLE, 
the mean square error should behave as 



V N {M, 9, 9 M le) = H{9)- 1 /N + o(l/N). 



IV. OPTIMAL INPUT STATE 



(0) 



It is now clear from Eq. © that one needs to choose 
the input state so that the QFI is "large." Suppose there 
is an input state that has a QFI that is larger than or 
equal to the QFI of any other input state. In that case 
one can unambiguously choose that state as the optimal 
one. Unfortunately, in our case there is not such a state. 



can be taken to be pure since it was shown in that the QFI 



Since the self-adjoint matrices T±, . . . , Td-i commute 
with one another, there is a basis where all of them 
are diagonal {|1), . . . , \d)} (this basis is considered to be 
known). From this point, all calculations will be made in 
this basis. 

The input state is po — \^o)(^o\, and |^o) can be 
expanded as |^o) = Ylu Rkl\kl)- The partial trace is 
then tvs |^o) (^o| = RR, and since RB) is self-adjoint 
and has trace 1, it can be written as 



RW 



1 
d 



where, in general, the last sum includes all generators of 
the Lie algebra su(d). Then the QFI is 



(7) 



H mn = 4[(tf |(T m r„<8>l)|*o> 

- (*o|(T m ® l)|* )(*o|(7 1 „ ® l)|*o>] 
= 4 [tr(RR^T m T n ) - (tr RR)T m ) (tr RR?T n )] 
= 4 [tr(RR)T m T n ) - t m t n ] 
and its trace is 

(d-l \ d-l 
Rrf T ™ - E 
m— 1 / m— 1 

The commuting T's can be written as T m — 
Sfc=i c mk\k){k\. The tracelessness condition implies 
Sfc=i c mk — 0; while the orthonormality condition (J5J 
implies Y^k=i c mkC n k — fimn- These two together lead to 



d-l ^ 
^ CmkCml — Ski ~ ~j 



and then to 

d-l 



is convex. 



rn — 1 k—1 m— 1 
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Substituting this in the equation for the trace, one gets 



TtH 



1 



and one immediately sees that tr H is maximal if and 
only if t m = 0, to = 1, . . . , d— 1. Of course the rest of the 
t's can be anything (as long as remains positive). 
Concluding, any pure input state po satisfying 

trp (T m <g> 1) = 0, V to = 1, . . . ,d - 1, 

achieves a QFI H mn = 4 S mn /d which has the maximum 
possible trace among all QFIs. In particular, maximally 
entangled states satisfy this condition. 

In the full SU(d) model one obtains a similar result, 
and the maximum is attained at and only at the max- 
imally entangled state; this will be shown in appendix 

HI 



B. MPEU 

It will be shown here that for every entangled input 
state |* > = Y,ki R ki\kl) eC^C d and every value of 
the [d — l)-dimensional parameter 9 there is an input 
state \ipo) G C d that achieves the same QFI everywhere. 
This means that in this model, unlike in the full SU(rf) 
model, there is no improvement in the accuracy of the 
estimation by using entangled inputs. 

The model is now 



p{6) = UepoUl 



where 



/ d-i \ 
U e = cxp i 22 m T m 
\ m=l / 

and po = | -00 ) {^o I is a pure state density matrix. Of 
course this model is also quasiclassical and by a calcula- 
tion identical to the one made in the C d <8> C d case one 
gets that the QFI is 

\ipo) - (V'o|T 1 m |V'o)(V'o|'7 1 n |V'o)] ■ 
It is not difficult to see that any input state of the form 



d 

E 

k=l 



(k\RRt\k) e^ k \k) 



where the </>'s are arbitrary phases, achieves 

ff m „ = 4p |(T m T n (8l)|*o) 

- {* \{T m <g> l)|*o)<*o|(T n ® i)| * )] , 

the QFI when the input is the entangled state \^>q) 
In particular, any input state of the form 

d 



1^0 



-y 



■\k) 



(8) 



k=l 



achieves the maximum trace of the QFI. A state of this 
form (with all the </>'s set to zero) was used in Ref. 0. 



V. OPTIMAL MEASUREMENT 
A. MPEU 

This model (and actually also the MPEE) has the 
property that I(6,UgMUg) = 1(0, M), which is easy to 

prove; the only necessary ingredient is that Ugdg k Ue — 
d$ h Ug\g—o- Therefore, if one has an optimal measure- 
ment M at 9 = 0, then the measurement UgMUg will be 
optimal at 9. 

One can now use the recipe given by Eq. @ to find an 
optimal measurement at 9 = 0. The optimal input state 3 
given by Eq. JSJ is used. At the origin one has 



2i 



^E c «^ e *l fc >' 

71 



fc=l 

From these vectors one can form an orthonormal set 



\m„ 



= i E Cnk e ^ fc n = l,...,d-l, 



k=l 



|m d ) = ^=E e ^l A; ) 
^ d k=i 

and with the choice Oki = Ski — 2/d one gets the set 
of orthonormalized states onto which the measurement 
elements will project, 



\b k ) = \m k ) - -^2\mi), k=l, 



(9) 



i=i 



The optimal measurement at 9 = has elements Mk = 
\bk)(bk \ and since the above vectors form an orthonormal 
basis of C d , they satisfy ^fe=i = 1- The optimal 
measurement at 9 has elements Ue\bk){bk\Ul . 

Note that the above choice of the orthogonal matrix 
o works for d > 3; for d = 2 another matrix must be 
chosen. The d = 2 case will be treated in the following 
example. 

Example V.l (d = 2) The orthonormal set formed 
from the input state and \l) is 



\mi) = -j=(\Q) |1» 
| m2 ) = -L(|0) + |l». 



(10) 



the one that achieves maximum trace of the QFI 
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These vectors are then rotated to obtain 



VI. ASYMPTOTIC FIDELITY 



1 61) = cos?7|mi) - sin77|m 2 ) = -r=(e in |0) - e "'|1)) 

f (11) 
|6 2 > = sin?7|mi) + cos77|m 2 ) = —={e lv \Q) + e~"'|l)), 

v2 

where rj must satisfy sin 77 7^ and cos 7/ 7^ 0. T7ie mea- 
surement elements are \b\)(bx\ an d I^K^I- The FI of 
this measurement at the origin is equal to the QFI ( equal 
to 2) as expected. Furthermore, this equality happens to 
hold everywhere and not only at the origin. This feature 
is very useful in practice, it means that the optimal mea- 
surement does not depend on the actual (unknown) value 
of the parameter and therefore an adaptive scheme is not 
necessary. Whether it is possible to find measurements 
with this characteristic for d > 2 is an open problem. 



B. MPEE 

An optimal measurement in this case can also be de- 
rived using the recipe given by Eq. |J5J. In general, such 
a measurement is a joint measurement on the two output 
systems. One could ask whether it is possible to achieve 
the bound with local measurements and classical com- 
munication. In what follows, it will be shown that this 
is indeed possible. 

The input state is taken to be the maximally entangled 
state J2k=i \kk}/Vd. Alice measures the system coming 
out of U and Bob measures the other one. The strategy 
is the following. Bob performs the von Neumann mea- 
surement B k = \wk)(wk\ with J2t Bk — 1 on his system 
where 



1 

7^ 



exp 



{=1 



2irkl 



i |0, k = l, 



He obtains outcome k with probability l/rf; he then 
phones Alice and tells her the outcome of his measure- 
ment. The net result of this is that Bob prepares the 
state 



1 

— = exp 

Vd^ 



2irkl 



10 



at the input of U . It is easy to recognize this state as 
one of the optimal states in the C d case. Now Alice can 
perform the measurement A^i with A^i = 1, which 
is the optimal measurement described above for the C d 
case, and where the arbitrary phases are now fixed to 
<f>i = 2-Kkjd. It is not difficult to check that this mea- 
surement indeed achieves equality in the QCRB at 9 = 0. 
The measurement on C d ®C d is then ^2 kl A M ®B k — l<g)l 

and the measurement ^2 k[ UsAkiUg (g> Bk = 1 ® 1 is op- 
timal at 9. 

If this is applied to the d — 2 case and one uses the 
^-independent optimal measurement derived for MPEU, 
a ^-independent optimal measurement is obtained for 
MPEE. 



From the results obtained previously, one can infer 
the asymptotic behavior of the average fidelity. Indeed, 
in Ref. [jg, it was established that the fidelity between 
nearby states is given by 4 



T{9, 9 + 89) = ( tr sj '^(9) p(9 + 59) ^(9) 



^^59 a S9 p + o(S9 2 ), 

a, /3=1 



(12) 



where p is the number of parameters and H is the QFI. 
In the case studied here, this fidelity would be between 
output states. 

Denote by 9% the guess for 9 if the outcome of the 
measurement was £, then the fidelity, averaged over all 
possible outcomes, is 



= X>b(0)M € ] no A) 

tiH(9)V(M,9,9) 



(13) 



1 



o(S9 2 ), 



where V(M, 9, 9) a0 =Y,^W) M t\ {0 ta -9 a ){9 tp -9 p ) 
is the mean square error of measurement M and estima- 
tor 9. 

Therefore, using Eq. |(BJ| one gets 



iV 



lim N[l - F(6»,6») JV )] = 



Tr H(9)H(9Y 



d-l 



(14) 



This result can be compared with the optimal fidelities 
obtained in Ref. 6] . These optimal fidelites were also av- 
eraged with respect to a uniform prior distribution on 9. 
However since the result obtained here does not depend 
on 9, its average with respect to any prior will be itself. 
The comparison is shown in Fig. [3J 

One can easily see that for large N the optimal fidelity 
of Ref. agrees with the result obtained here. Actually, 
this can also be proved analytically but the proof will not 
be shown here. 



VII. CONCLUSIONS 

Two models have been compared, the model of esti- 
mating commuting unitaries with and without the use of 
entangled inputs (MPEE and MPEU, respectively). 

It has been shown that the quantum Cramer-Rao 
bound is achievable in both MPEE and MPEU. It has 
also been shown that any quantum Fisher information 



4 Actually, this fidelity is the square of the fidelity used in Ref. |^ . 
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FIG. 3: The points are N(l - F N ) as a function of N for 
d = 2, . . . , 5, where F N is the optimal fidelity obtained in 
Ref. |g. The continuous lines are at the value (d — l)/4 for 
d = 2, . . . , 5, (0.25, 0.5, 0.75, and 1, respectively.) 



matrix that can be attained in MPEE can also be 
achieved in MPEU. These two facts imply that an entan- 
gled input state is unnecessary. A condition for attaining 
maximal trace of the QFI has been derived. 

In the MPEE it has also been shown that there is a sep- 
arable measurement that achieves equality in the QCRB. 

In the d = 2 case, measurements that are optimal ev- 
erywhere have been found in both MPEU and MPEE. 
This is a useful feature in practice since this means that 
an adaptive scheme would not be necessary in this case. 
However, it is unclear whether it is possible to find mea- 
surements with this characteristic in general. It is also 
an open question whether entanglement could prove itself 
useful in this respect (for d > 2). 

These facts show that entanglement is, at best, not as 
useful in estimating commuting unitaries as in the esti- 
mation of a completely unknown unitary. 
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APPENDIX A: GENERALIZATION OF THE 
RESULT IN SECTION [TV] TO THE FULL MODEL 



In this appendix, a result similar to that of Sec. IIVI will 
be proved in the model that includes the whole of SU(d) 
and not only a commuting subgroup, i.e., the model con- 
sidered in Ref. 

Denote by H Po (9) the QFI at 9 if the input state is poj 
and by H the QFI when the input state is a maximally 
entangled state; in what follows the dependence on 9 will 



be omitted. Then the following inequality holds for any 
input state po : 



Tr(H- 1 H p0 ) < d 2 - 1, 



(Al) 



and equality is attained if and only if po is a maximally 
entangled state. This will be proved in what follows. 

Notice that this trace is parametrization invariant, and 
not just Tr H p as in Sect. IIVI in that case H was propor- 
tional to the identity so there is no contradiction. 

The model is again described by Eq. Q but now U is 

U e = V e ® 1 



where Vg = cxp [iJ2 a =i 9 a T a J. As before, the T's arc 

traceless self-adjoint matrices chosen so that tr(T a Tp) = 
5 a g. The input state po is chosen to be pure because 
of the convexity of the QFI 8]. The SLDs are X a (9) = 
2p, a {9), where p, a (9) means the partial derivative of p(9) 
with respect to 9 a . The matrix elements of H Po are 



ttPo 



Retr[pA Q A^] = 4Retr[/jp iQ p i( g]. 



Since po is pure it can be written as po — |V'o)('0o| an d 
l^'o) = Xlfcz Rki\kl)- H^°g(9) can then be calculated to be 

fl$ = 4Re [tr^tyt^) 

+ triRR^V^V^triRR^Vp)] . 

Denote S a — —iV^Va; then 

= 4 Re [tr(RR*S a Sp) - tr(RR*S a ) tr(RR^Sp)] . 

Note that S a £ su(d). Substituting RR^ = t/d in the 
expression for H Po ; one gets 

4 

H a (i — — ti(S a Sp). 
The matrices Si, . . . , S^-i can be orthonormalized, 



tr 



The operator 



is a Casimir operator and therefore proportional to the 
identity. The proportionality factor can be found by tak- 
ing the trace, and finally one gets 



d 2 -l 



t. 
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The wanted trace is 
TriH^H" ) = d 2 -l-^2 tr(RR?S a ) tr(RR^Sp). 

a/3 

This quantity is always less than or equal to d 2 — 1 
and, furthermore, this value is attained if and only if 



tr(RR^ S a ) — for all a = 1, . . . , d 2 — 1, which implies 
that RRJ = i.e., po is maximally entangled. In par- 
ticular, this implies that there is no input state po for 
which H Po > H. 
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